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ABSTRACT 

We investigate the thermal stability of optically thin, two-temperature, ra¬ 
diative cooling-dominated accretion disks. Our linear analysis shows that the 
disk is thermally unstable without magnetic helds, which agrees with previous 
stability analysis on the Shapiro-Lightman-Eardley disk. By taking into account 
the effects of magnetic helds, however, we hnd that the disk can be or partly be 
thermally stable. Our results may be helpful to understand the outhows in opti¬ 
cally thin hows. Moreover, such radiative cooling-dominated disks may provide 
a new explanation of the diherent behaviors between black hole and neutron star 
X-ray binaries on the radio/X-ray correlation. 

Subject headings: accretion, accretion disks - black hole physics - instabilities - 
magnetic helds - stars: winds, outhows 


1. Introduction 

Accretion of rotating matter onto a compact object can provide a large amount of 
released energy and is therefore believed to be the source of cataclysmic variables. X-ray 
binaries, and active galactic nuclei (AGN). The most famous accretion model is the geo¬ 
metrically thin and optically thick disk (SSD) introduced by Shakura & Sunyaev (1973). 
Thermal instability was found in the inner region of SSD (Piran 1978) and such an issue has 
been widely studied by linear analysis and numerical simulations (Hirose et al. 2009; Lin et 
ah 2011, 2012; Xue et al. 2011; Jiang et al. 2013; Zhu & Narayan 2013). For stellar-mass 
black holes, the temperature in such a standard disk is in the range 10^ — 10^ K, which is 
quite low relative to the virial temperature. In order to explain the hard X-ray spectra of 
black hole X-ray binaries such as Gyg X-1, Shapiro, Lightman & Eardley (1976) introduced 


an optically thin, two-temperature disk model (hereafter SLE disk) with Te ~ lO^K. Unfor¬ 
tunately, shortly after the SLE disk was introduced, it was found that the disk is thermally 
unstable (Pringle 1976; Piran 1978). Both of the above two models are radiative cooling 
dominated, i.e., the advective energy is negligible. 

On the other hand, the energy released through dissipation may be trapped within the 
accreted gas and then transported in the radial direction towards the central object or stored 
in the flow as entropy. Such accretion flows can be divided into two types, namely the slim 
disk and the advection-dominated accretion flow (ADAF). The slim disk was introduced by 
Abramowicz et ah (1988), which is optically thick with extremely high mass accretion rates. 
Recently, simulations on the slim disk have revealed that strong outflows may occur (e.g., 
Ohsuga et al. 2005; Ohsuga & Mineshige 2011; Yang et ah 2014). Moreover, the radiative 
efficiency is not low and the luminosity can be far beyond the Eddington one (Jiang et al. 
2014). The ADAF model was proposed by Narayan & Yi (1994), which is optically thin 
with low accretion rates. The temperature of ADAF is close to the virial one, which is 
signihcantly higher than that in the SLE disk due to energy advection. Simulations on the 
ADAF have been studied by many previous works (Narayan et al. 2012; Yuan et al. 2012a,b; 
for a review see Yuan & Narayan 2014). Different from the stability of the SLE disk, both 
the slim disk and the ADAF are thermally stable due to the dominant energy advection 
(Abramowicz et al. 1995). A recent work (Gu 2014), however, showed that the advective 
cooling cannot balance the viscous heating, and therefore outflows ought to be inevitable in 
such an accretion system. 

The stability properties of accretion disks are of importance because a global violently 
unstable disk may not exist in nature and some instabilities restricted in a certain region of 
the disk may contribute to the observed light variations in many systems. Since the classic 
SLE disk was proved to be thermally unstable, such a type of accretion disk is unlikely to be 
realized in nature. It is generally believed that magnetic helds have fundamental influence on 
the physics of accretion disks. For example, magnetorotational instability (MRI) is known as 
a generator of turbulence, through which the angular momentum can be transferred outwards 
(Balbus & Hawley 1991). Oda et al. (2009, 2010) presented a new thermal equilibrium 
solution for optically thin disks by incorporating toroidal magnetic helds. They argued that 
an optically thin, magnetic pressure-dominated accretion disk (low-/? disk) will be thermally 
stable. Zheng et al. (2011) revisited the thermal stability of standard thin disks by including 
the role of toroidal magnetic helds. Their calculation is based on the assumption of the 
perturbation relation 5B^p/B^p = —'ySH/H, where 7 is positive, i.e., the magnetic held will 
become weaker with increasing height (or temperature). This assumption is supported by 
the MHD simulation of Machida et al. (2006) for a hot accretion how, where they found 
that the magnetic held becomes stronger when the disk shrinks vertically. In addition, the 
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large-scale magnetic fields in accretion disks have been widely investigated and the strength 
may be underestimated particularly for the region near the black hole (Cao 2011). 

In this paper, we will analyze the thermal stability of optically thin disks by including 
the role of magnetic fields. The effects of such magnetic helds can be taken into account 
by modifying the pressure as the sum of the gas pressure and the magnetic pressure: p = 
Pgas + Pmag, Pmag ^ ^ (Narayau & Yi 1995; Yamasaki 1997). We have to consider 

two-temperature plasma in such disks because the electron temperature is expected to be 
signihcantly lower than the ion temperature in such a low density, high temperature region. 
In addition, the radiative cooling process is assumed to be the thermal bremsstrahlung. The 
paper is organized as follows. The basic equations are presented in Section 2. In Section 3, 
the thermal stability is investigated by linear analysis. In Section 4, the thermal equilibrium 
curve is obtained by numerical calculation. Conclusions and discussion are made in Section 5. 


2. Equations 

We consider a steady state, optically thin, two-temperature black hole accretion flows 
incorporating magnetic helds. We assume the angular velocity is Keplerian, i.e., hi = IIk- 
The disk structure is described by the following equations. The continuity equation takes 
the form: 


M = —AnRHpvR , 


( 1 ) 


where R is the cylindrical radius, p the mass density of the accreted gas, H the half-thickness 
of the how, vr the radial velocity, which is dehned to be negative when the how is inward, 
and M the mass accretion rate. The equation of vertical hydrostatic equilibrium can be 
written as 


H 


1 



( 2 ) 


where p is the pressure. The Keplerian angular velocity fix takes the form = GM/[R — 
Rg)^R, where the gravitational potential of the central black hole is assumed to be *h(7?) = 
—GM/[R — Rg), which was introduced by Paczyhski & Wiita (1980), with M being the black 
hole mass and Rg the gravitational radius, Rg = 2GM/c^. From the angular momentum 
equation we can obtain the expression of the radial velocity (e.g., Gu et al. 2006; Liu et al. 
2007), 


vR/^ dVt 

VtR^ -jdR ’ 


Vr = 


( 3 ) 
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where v is the kinetic viscosity coefficient, and j represents the specihc angular momentum 
per unit mass accreted by the black hole. 

Since the radiation pressure is negligible in optically thin disks, the total pressure p is 
expressed as 


V = Pgas + Pmag , (4) 

where Pgas = (1 —/dmag)p is the gas pressure, pmag = /dmagP the magnetic pressure. The gas is 
assumed to consist of protons and electrons and therefore the gas pressure Pgas is written as 

Pgas = -^p(Ti+Te) , (5) 

prUp 

where ks is the Boltzmann constant, p, the mean molecular weight {p = 1), nip the proton 
mass, and Tj, Tg the ion and electron temperature, respectively. 

To construct the energy equations, we introduce the following assumptions: i) the dis¬ 
sipation energy goes into the ions; ii) the energy is transferred from the ions to the electrons 
through the Coulomb coupling; iii) the electrons are cooled by the bremsstrahlung process. 
Based on these assumptions, the energy equations of the ions and electrons can be respec¬ 
tively written as 


= ^ie , ( 6 ) 

^ie — Qe ) ('^) 


where is the viscous heating rate per unit volume, Aig the energy transfer rate from the 
ions to the electrons per unit volume, q~ the radiative cooling rate of electrons per unit 
volume. The viscous heating rate can be written as 


We adopt the standard a prescription in this paper, i.e., u = aCsH, where Cg = {p/pY^"^ and 
a is a constant parameter. Ajg is given by Stepney & Guilbert (1983): 


3 mg 2 n kTi-kT^ 

^ (In A) ^ 

2 mi A2(l/Bi)A^2(l/0 

2(0i -|- 0e)^ + 1 f Qi + Q, 


0 i -|- 0 g 


-K, 


0i0e 


+ 2A’o 


0i -|- 0e 
0i0e 


(9) 


where is the Thomson scattering cross section and In A is the Coulomb logarithm (roughly 
InA ~ 20). Kn are modihed Bessel function of the second kind of the order n, respectively. 
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The quantities ©i and ©e are defined as ©i = k-BTi/m^c^ and ©e = k-BTe/rrieC^, where rrip 
and nie are the mass of the proton and electron, respectively. For simplicity, we use the 
following formula which uses no special functions (Kato et ah 2008): 


3 p/CB(Ti-Te) 

Aie — -r'E- , 

2 ITbp 

with the electron-ion coupling being z/e = 2.4 x 10^^(ln A)pTe”^'^^. 


( 10 ) 


Following Narayan & Yi (1995), the bremsstrahlung cooling rate per unit volume is 

C = = Cr,ei + Cr,ee = (TTCafTTleC^ [F;i(©e) + F;e(©e)] , (H) 


where the subscripts “ei” and “ee” denote the electron-ion and electron-electron bremsstrahlung 
cooling rates, ctf is fine-structure constant, and the function Fei(©e) and Fee(©e) have the 
approximate form: 


r 9©. 


27r 


ln(2r/©e + 0.48) + 


Fei(©e) = < 


i"ee(©e) = < 


2©e 

9©e 

71 

5 


1/2 


67J-3/2 


(©e > 1) 


[l + 1.781©^-^^] (©e < 1) 

[ln(2p©e) + 1.28] (©e > 1) 

(44 - 37r2)©px 


[ (1 + l.l©e + ©^ - 1.25©e''") (©e < 1) 

where t] = exp(— 7 e) and 7 e ~ 0.5772 is Euler’s number. 


( 12 ) 


3. Thermal stability analysis 


The thermal instability criterion of accretion disks can be expressed as (e.g., Frank et 
al. 1992) 


/(91n(5 \ ^ /(91n(5’''\ 


(13) 


where , Q~ and T are respectively the viscous heating rate, radiative cooling rate, and 
temperature. S is the surface density defined as S = 2pH. Such a criterion is identical to 
that in Piran (1978) if we replace T by H. 
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Since the disk is two-temperature, the thermal stability analysis should be made for ion 
and electron temperature modulations, respectively. The cooling timescale of the electron is 
quite short and therefore the thermal stability of the electron is priority. The surface density 
S is assumed to be unchanged during a thermal timescale, so we have 

d In p d In id 


dluTp 


dluTp 


(14) 


With Equations (7), (10), (11) and (14), and the assumption that Ti keeps unchanged during 
the perturbation of Tg, we can obtain 


/ (91ng+ 


/ (91ng^ 

\dlnT, 


< 0 , 


(15) 


e / E 


where is equal to Ajg. The above relationship means that the disk is always stable against 
electron temperature perturbations. Such a result is in agreement with Kato et al. (2008). 


If electrons are thermally stable, then we have 

d(lng+) d(lng; 


dluTj dluTj 


(16) 


With Equation (7) we can obtain the ratio of T, to Tg as a function of Tg, i.e., a function of 

0e, 


Equation (2) gives 


Y = sm ■ 


dlnp dlnp 

dluTi dlnTj 


(17) 


(18) 


By using Equations (4), (5) and (18), together with the relationship Pmag oc p^/^ mentioned 
in Section 1, we can derive 


dlnp _ 3(1 - /^mag) dln(Ti Tg) 

d In Ti 6 -|- /djnag d In T, 

then with Equations (16), (17) and (19) we can obtain 

dTp 


dTi 

Equations (6), (8) and (10) can provide 


— h(/3niag) ©e) • 


dln(gj^ —gj) ^dlnp dln(Ti —Tg) ^ 3dlnTg 
= “d-r; — — — h —- 


dluTj 


dluTj 


dluTi 


2 d In Tj 


(19) 


( 20 ) 


( 21 ) 



7 


Substituting Equations (17), (19) and (20) into Equation (21), we finally get 

= /(/9mag,0e) , (22) 

S 

then the thermal instability criterion can be expressed as /(/Smag, ©e) > 0. Detailed expres¬ 
sions of g{Qe), h{/3ma.g, ©e), and /(/Smag, ©e) are presented in the Appendix. 

Figure 1 shows the variation of the function /(©e) for a certain given value of /Smag- It 
is seen that for /3ma.g = 0, i.e., no magnetic field, /(©e) is positive for any 0e, which means 
that the disk will be thermally unstable. This result agrees with previous stability analysis 
on the SEE disk. Actually, for /Smag = 0, the equations will reduce to that for SEE disks. 
The interesting result is that for /Smag = 0.7, /(©e) is negative for any ©g, which implies that 
the disk will be thermally stable. For the typical equipartition case /3ma.g = 0.5, the figure 
shows that there exists a thermally stable region corresponding to 0.03 ^ ©e ^ 1, whereas 
for ©e < 0.03 or ©e > 1 the disk will be thermally unstable. In other words, the disk can be 
or partly be thermally unstable owing to the effects of magnetic fields. 

The physical interpretation of thermally stable disks for moderate strength of magnetic 
fields may be as follows. For a positive perturbation of the ion temperature Ti, the gas 
pressure Pg^s will increase and therefore the disk thickness H will also increase. Consequently, 
the strength of magnetic fields B and the magnetic pressure Praa.g will decrease due to the 
increasing volume. Thus, the total pressure {pgas + Pmag) as well as the viscous stress will 
not increase as fast as in the disk without magnetic fields. Therefore, the viscous heating 
rate may increase more slowly compared with the case without magnetic fields. In this 
spirit, it is reasonable that the disk can become thermally stable when the role of magnetic 
fields is taken into account. 


(91n(g+ - q. ) 
(91nTi 


4. Thermal equilibrium solutions 

In this section, we will show the numerical calculation of the local thermal equilibrium 
solutions with Equations (l)-(7), where M = IOMq, j = 1.8ci?g, a = 0.1, and /Smag = 0.5. 
Figure 2 shows the radial variations of the ion and electron temperature for m = 0.001, 
0.01, and 0.1, where m = M/MEdd {M^dd = 16EEdd/c^)- It is seen that the ion temperature 
is significantly higher than the electron temperature in the inner region of the disk. We 
would like to point out that the temperature in the real case may be lower than the current 
solutions since only the bremsstrahlung cooling process is taken into consideration. In the 
classic ADAF model where the radiative cooling is negligible compared with the viscous 



heating, there exists a conservation relationship (e.g., Molteni et ah 2001), 


Vr + 


7 P 


7 - Ip 


+ - j) - 


GM 
R — Ra 


= A, 


( 23 ) 


where A is a constant and should be quite small compared with the terms on the left hand 
side such as GM/{R — Rg), particularly for small radii. Based on the above equation, we 
can derive the theoretical ion temperature of ADAF by 


5 /ceTj 

2 /imp 


GM 
R — Rg 


+ -^^R^ — jfl , 


(24) 


where 7 = 5/3 and A = 0 are adopted, and the term n^/2 is ignored. The theoretical ion 
temperature is also plotted in Figure 2 by the navy blue line for a comparison. It is seen 
that the ion temperature of ADAF is significantly higher than the cooling-dominated disks 
even for m = 0.1, which indicates that the energy advection is negligible and our thermal 
equilibrium solutions are self-consistent. 


Figure 3 shows the stable and unstable regions in the m — R diagram for /dmag = 0.5. 
The region with /(0e) < 0 corresponds to the thermal stable solutions. For a given m, the 
hgure clearly shows the thermally stable range in the disk. It is seen that, for a typical 
accretion rate m = 0.01, the disk is thermally stable from the inner boundary to IQ^Rg. 
Such a stable solution may provide a second possibility for the optically thin accretion flows. 
For higher accretion rates such as m > 0.1 the hgure shows an unstable inner region, which 
implies that outhows are likely to be inevitable and the rate accreted by the black hole is less 
than 0.1. On the other hand, for low accretion rates such as m < 0.001, the disk only has an 
inner stable part < 100i?g. A thermally stable disk for such low accretion rates may require 
stronger magnetic helds, i.e., /Smag 0.7. Another possibility for the outer part of the disk 
to be stable may be related to the self-gravity, as argued by Bertin & Lodato (2001). 


5. Conclusions and Discussion 

In this work, we have investigated the thermal stability of magnetized, optically thin, 
two-temperature, radiative cooling-dominated accretion disks by linear analysis. We have 
derived a general criterion for such an instability (Equation 22). We have found that the disk 
is thermally unstable without magnetic helds, which agrees with previous stability analysis 
on the SLE disk. On the contrary, for adequately strong magnetic helds with /Smag ~ 0.7, 
the whole disk will be thermally stable. For the typical equipartition case, /Smag = 0.5, the 
disk has a wide stable region particularly for m around 0.01. Our results may be helpful 
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to understand the mechanism of outflows in optically thin accretion flows, which may be 
triggered by the thermal instability. 

The well-known ADAF model is a classic model, and may be known as the unique stable 
model for optically thin accretion flows. Since the black hole has a horizon rather than a 
solid surface, it is easy to understand that the flow will be transonic and the internal energy 
will eventually be absorbed by the hole. On the contrary, for a neutron star accretion system, 
since the neutron star has a real surface, it remains uncertain whether the ADAF can still 
be stable and work well in such a system, particularly for the inner region. The present 
work provides a second possibility for the optically thin accretion flows, i.e., the magnetized, 
radiative cooling-dominated disk. In our opinion, since the radiative cooling can balance the 
viscous heating locally in such a model, the energy advection will be negligible. Thus, such 
a model may work well for both black hole and neutron star systems. 

The observations have shown that the radio and X-ray emission during the hard state 
are strongly correlated for black hole and neutron star X-ray binaries. Taking the form 
of non-linear luminosity correlation, there exists the relationship Lr oc where Lr is 
the radio luminosity and Lx is the X-ray luminosity. For neutron star X-ray binaries, the 
correlation index is 6 ~ 1.4 (Migliari & Fender 2006). On the contrary, for black hole X-ray 
binaries, the correlation is quite complex. At hrst the correlation with b ~ 0.7 was obtained 
based on the data of different sources (e.g., Hannikainen et al. 1998, Corbel et al. 2003, 
Gallo et al. 2003), Later, however, another relation similar to that of neutron stars, i.e., 
b ~ 1.4, was found for specihc sources (e.g., Coriat et al. 2011, Jonker et al. 2012, Ratti 
et al. 2012). As discussed by Coriat et al. (2011), the steep correlation may indicate that 
the disk is in the radiatively efficient regime. In our point of view, the above observational 
results can be easily understood based on our new model. Due to the existence of surface, a 
neutron star system may prefer to the magnetized, radiative cooling-dominated disk which 
corresponds to 6 ~ 1.4, whereas a black hole system may have two choices, i.e., either the 
ADAF or the cooling-dominated one, which corresponds to & ~ 0.7 and ~ 1.4, respectively. 
In addition, we would point out that, an alternative cooling-dominated model for optically 
thin flows is the luminous hot accretion flow (LHAF) introduced by Yuan (2001), which is, 
however, likely to be thermally unstable. Concerning the stability, the magnetized, radiative 
cooling-dominated disk studied in the present work is more likely to exist than the other 
radiatively efficient disks. 
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A. Expressions of thermal stability criterion 

In this Appendix, we will derive the detailed expressions of g{Qe), h(/9jnag, ©e), and 
/(/dmag, ©e), which are mentioned in Section 3. 


fi'(©e) — 1 + 


[-Fei(©e) + -Fee(©e)] 


3.6 X 10^^ (mp + In A 


(Al) 


We dehne Xi, X 2 , Fi, 1^2, and k as follows. 


Xi = ©e 


ln(277©e + 0.48) + - 


+ 2©e [ln(2r/©e) + 1.^ 


X 2 = 8(2©e)^A(i + 1.78101-34) + (44 _ 37r2)0p(l + 1.10^ + ©e - 1.25©^A) ^ 

o 

Fj = [©e ln(2r/©e + 0.48) + 2r/©2/(2r/©e + 0.48) + 6©e + 2©e ln(2r/©e)] /Xi , 

Y 2 = [5.66©^-® + 37.1©^-®^ + 36©^-^ + 66©^'^ + 83©^'^ - 120©^] /X 2 , 

_ 3(1 - /jn,ag) 

6 ~\~ /3mag 

We can derive the following expression; 

l/{gm - 1) - 2k/{g{Q,) + 1) 


^(/3mag5 ©e) A 


1.5 + Yi + 1/((/(©e) ~ 1) + 2 / 1 ;/(5'(©e) + 1) 

l/(^(©e) - 1) - 2fc/(^(©e) + 1) 

1. 1.5 + 12 + 1/(5'(©e) ~ 1) + 2 / 1 ;/(^'(©e) + 1) 


(©e > 1) 
(©e < 1) 


(A2) 


Finally, we obtain 


/(/^mag,©e) 


3fc^(©e)(l + h(/jmag, ©e)) 

(fi'(©e) + 1) 


+ 1.55((©e)/i(/3mag,©e) 


^(©e)(l - h(/3mag,©e)) 

{gm - 1) 

(A3) 


which is the function introduced in Equation (22). 
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Fig. 2.— Radial variations of the ion temperatnre R (solid) and the electron temperatnre 
Te (dashed) for m = 0.001 (blue), 0.01 (green), and 0.1 (red), where /Smag = 0.5. The navy- 
blue solid line represents the theoretical ion temperature of ADAF (Equation 24). 













